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1 Introduction 

The Euclidean path integral formulation of quantum gravity is arguably the most promising 
of all non-stringy approaches to quantum cosmology. The canonical quantisation approach 
translated by the Wheeler-DeWitt equation does not possess the flexibility of the path inte- 
gral formulation as it cannot account for topology changing processes. The Euclidean path 
integral formulation with its sum over geometries naturally allows these processes to occur 
as it considers a sum over different topologies. As important as how to obtain the dynamics 
of the quantum state of the Universe is to use the correct prescription for the boundary 
condition needed to obtain the wave function of the Universe. The most natural boundary 
condition, the Hartle-Hawking no boundary proposal, |ll] , requires that the sum over geome- 
tries be restricted to those four-geometries which have no boundary except for the prescribed 
three-geometry where the arguments of the wave function are defined. This is specially well 
adapted to the path integral formalism as it can be easily implemented in a sum over topolo- 
gies context. However, even this proposal is plagued by at least two main problems. Firstly, 
the Euclidean gravitational action is not bounded from below which leads to the divergence 
of the Euclidean path integral. Secondly, there is no clear prescription for the correct inte- 
gration contour to use. In , Hartle proposes the use of the steepest descents contour in 
the space of complex metrics as the solution to both problems. Furthermore by choosing 
the steepest descents contour passing through the classical solutions of the theory, he made 
it very likely that the path integral be dominated by classical four-geometries, i.e., solutions 
of Einstein's equations and stationary points of the path integral, as desired for any wave 
function that is intended to represent our current Universe. Note that in this view, the fact 
that an integration solely over real-valued Euclidean geometries does not yield a convergent 
result for the path integral, is actually a good thing, for such a path integral would never 
predict the oscillatory behaviour in the late Universe that traditionally represents classical 
Lorentzian space-time. 

As it is clearly impossible to calculate the full wave function integrated over all possible 
metric degrees of freedom, it is usual to explore approximated models in which the infinite 
degrees of freedom are reduced to only a few. The loss of generality induced by this reduction 
is not as bad it might seem at first, if we remember the relative symmetry of the observed 
Universe. These models are generically called minisuperspace models, for they "exist" in 
small, few-dimensional, subspaces contained in the infinite-dimensional superspace which 



is the space of all three-metrics. Particularly useful among minisuperspace models are the 
the ones based on Regge calculus. Such simplicial minisuperspace models were introduced 
by Hartle [§]. In such models one typically takes the simplicial complex which models the 
topology of interest to be fixed and the square edge length assignments play the role of 
the metric degrees of freedom. The summation over edge lengths models the continuum 
integration over the metric tensor. This approach has several advantages. First by treating 
the four-geometry directly it is more adequate to deal with the Hartle-Hawking proposal, |1| , 
(with its four- dimensional nature), than the usual 3 + 1 ADM decomposition of space-time, 
where a careful study of how the four-geometry closes off at the beginning of the universe 
is essential. Second, by discretizing space-time the classical equations become algebraic 
which makes it easier to find classical solutions which are essential to the semiclassical 
approximation. Third, it offers the possibility of systematic improvement. 

In 0, Hartle studied the simplest simplicial minisuperspace model, the a^ triangulation 
of the three-sphere. He assumed isotropy which leads to all boundary edge lengths being 
equal. Later in [Q, Birmingham generalised his work to isotropic triangulations of Lens 
spaces L{p, 1). Another model was considered by Louko and Tuckey in |^ working in a 
minisuperspace anisotropic 2 + 1-dimensional cosmological model. In 0, Furihata considered 
a more general anisotropic triangulation of the three-sphere. Finally, in |^, we generalised 
the isotropic models to incorporate all simplicial four-geometries that are cones over closed 
connected simplicial three- manifolds, i.e., simplicial four-conifolds, and for the first time 
introduced a massive scalar field. In all these cases it was possible to find a steepest descent 
contour and to prove that the resulting wave function shows oscillatory behaviour for large 
spatial geometries. We now propose to study an extension of Furihata's model to incorporate 
a matter sector. We investigate the existence of a steepest descent contour and the properties 
of the resulting wave function. 

2 General Regge Formalism 

A convenient way of defining an n— simplex is to specify the coordinates of its (n + 1) vertices, 
a = [0, 1, 2, ..., n]. By specifying the squared values of the lengths of the edges [i,j], Sij, we 
fix the simplicial metric on the simplex. 

g^A^k) = ''' ^ ''' ~ ''' (2.1) 



where i,j = 1, 2, ..n. 

So if we triangulate a smooth manifold M endowed with a metric g^i, by a homeomorphic 
simplicial manifold Ai, the metric information is transferred to the simplicial metric of that 
simplicial complex 

9.A^) — g^Aisk}) = "'' ^ 'I' ~ "'' (2.2) 

In the continuum framework the sum over metrics is implemented through a functional 
integral over the metric components {g^u{x)}. In the simplicial framework the metric degrees 
of freedom are the squared edge lengths, and so the functional integral is replaced by a simple 
multiple integral over the values of the edge lengths. But not all edge lengths have equal 
standing. Only the ones associated with the interior of the simplicial complex get to be 
integrated over. The boundary edge lengths remain after the sum over metrics and become 
the arguments of the wavefunction of the universe. 

j Dg^,{x) — . j D{si] = Y[Jdfi{si) (2.3) 

In the simplicial framework the fact that the geometry of the complexes is completely 
fixed by the specification of the squared values of all edge lengths, means that all geometrical 
quantities, such as volumes and curvatures, can be expressed completely in terms of those 
edge lengths. Consequently the Regge action (the simplicial analogue of the Einstein action 
for GR) associated with a complex of known topology can be expressed exclusively in terms 
of those edge lengths. 

The Euclideanized Einstein action for a smooth 4— manifold M with boundary dM, and 
endowed with a 4— metric, g^^, and a scalar field $ with mass m, is 

I[M,g,^,<j>] = - I d'x^g^^-^ - j d'xVh^ + 
Jm IbTrG JdM SttG 

+ 11 rf^xV^(c>^0a'^0 + mV) 
2 Jm 

where K is the extrinsic curvature. 

Its simplicial analogue will be the Regge action for a combinatorial 4— manifold, A4, with 
squared edge lengths {sk}, and with a scalar field taking values {4>v} for each vertex v of 
M, [HI: 
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+ ^E^4(jK0? 

where: 

• (Tfc denotes a /c— simplex belonging to the set S^ of all /c— simplices in M.. 

• Oicl), is the deficit angle associated with the interior 2— simplex a^ = [ijk] 

9ial)=2n- E ^Ml^^^) (2-4) 

and 9fi{a\, (j/i) is the dihedral angle between the 3— simplices a^ = [ijkl] and a^ = 
[ijkm], of o"4 = [ijklm] that intersect at al. Its full expression is given by [Q. 

• ip{(J2) is the deficit angle associated with the boundary 2— simplex o"2: 

,p{4) = n- E ^<i(^2,^4) (2.5) 

• Vfc(crfc) for fc = 2,3,4 is the A;— volume associated with the fc— simplex, ak, and once 
again their explicit expressions in terms of the squared edge lengths are given by Q]. 

• V4(o"i) is the 4— volume in the simphcial complex JH, associated with the edge cxi, i.e., 
the volume of the space occupied by all points of Ai that are closer to ai than to any 
other edge of A4. The same holds for Vi{j) where j represents all vertices of A4. 

It is easy to see that both V^^ai) and V^i^j), can be expressed exclusively in terms of the 
edge lengths {sk}- All these expressions remain valid if we consider smooth conifolds and 
their combinatorial counterparts. 

If we are working in a simplicial minisuperspace model based on a simplicial four-manifold 
A^4, with boundary edges Sb, interior edges Sj, and in the presence of a scalar field taking 
values, 0fc, at the vertices k G A^4, then the wave function of the Universe is of the type: 



^OM, {sb}, {(j)b}] = f D{si}D{(j)i}e-'^^''^'^^'^''^'^'f'^^'^^'^^ (2.6) 

where 

• {si} are the squared lengths of the interior edges 

• {sb} are the squared lengths of the boundary edges 

• {(pi} are the values of the field at the interior vertices 

• {0;,} are the values of the field at the boundary vertices 

Although the functional integral over metrics has been written explicitly in terms of the 
edge lengths, this expression is still heuristic because we still need to specify the measure, 
and the integration contour to be used. 

3 The Model 

As in Q our simplicial complex can be seen as a cone, Ai"^ = * Ai^, over an anisotropic 
triangulation of the three-sphere, A4^ = 04(0, b). See figure 1. There are 10 boundary edges, 
five with squared length, a, and the other five with squared length b. There are also five 
interior edges connecting the five boundary vertices 1,2,3,4,5 to the interior vertex 0, and 
they all have squared length Sj. Note that although the model is anisotropic because there 
are two kinds of boundary edges, it is still homogeneous because all boundary vertices are 
connected to the same number and kind of edges. Furthermore, we consider the scalar field 
to take the same value, (pb at all boundary vertices, and (pi at the interior vertex 0. 
For simplicity we will use new variables. 
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The Euclidean action /, is then a function of ^, S, (3, (pb, 4>i, and the minisuperspace wave 
function is 

^[S,(3,(j)b] = f D^D0ie-^K'^''3'<^-<^^l (3.4) 

To calculate the explicit expression of the action we must compute the volumes and 
deficit angles associated with the various p — simplices in A^4. 

There is only one kind of 4 — simplex, eg., [01234], and its volume is 

1 /a 



V401234] = -i^-j ^-/52 + 3/3 -lv/4e(/5+l)- 1-/5-/52 (3.5) 

There is also only one kind of boundary tetrahedron, eg., [1234], and its volume is 



^3[1234] 



2 /a 



(0'(/3+l)(-/9^ + 3/3-1) (3.6) 



The existence of two different boundary edges means that there are two kinds of interior 
tetrahedra 

^3a[0123] = (i) {^^Y'^2P{{A-P)^-1} (3.7) 

V^3'i.[0124] = (i) g)'^V2{(4/5-l)e-/3^} (3.8) 

There are also two kinds of interior triangles, whose areas are 



a 



^2A[012] = -v/e-l/4 (3.9) 



V^2B[013] = ^(e-f)'^' (3.10) 



Finally there are also two kinds of boundary triangles 



V'24[1241 = -\//3 - 1/4 (3,11) 

ViBm = ^{l-!^f (3.12) 

The 4 — volumes associated with each boundary vertex b, and with the interior vertex 
are 



VM = -\/4[01234], VM = \/4f01234] 
5 

The 4 — volume associated with each interior edge [0, b] is 



(3.13) 



^4(0, 6] = -l^4[01234] 
5 

The deficit angles associated with each kind of triangle are 



(3.14) 



^[012] = 27r — arccos hi — 2 arccos/i2 
^[013] = 27r — arccos h^ — 2 arccos /i4 
= TT — 2 arccos h^ 



where the dihedral angles are 
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h: 



2p{A^ - 1)(4 - /?) v^(e - i,A){i - i,B) 
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where 
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/3^ 



4-/3 """ 4/3-1 

are the values of ^ for which the interior tetrahedra become degenerate. 



(3.25) 



The Regge minisuperspace action then becomes 



/[e,5,/3,0„0,] = -(-^)|-V^4/3 - P\ti - 2arccos/i6) + -^A(3 - l(7r - 2arccos/i 



+ bJ^ — l/4(27r — arccos hi — 2 arccos h 



2 



+ hJ j3\l^ — -ri^Ti — arccos /i3 — 2 arccos /i4) 



- ^y^nv-/?^ + 3/3 - 1^^^(0. - c^,ff] 



where 



(3.26) 

^' - 4/5 + 4 ^^-^^^ 

is the the value of .^ for which the 4 — simplices become degenerate. 

Since the wave function is to be obtained as an integral over C, and (pi, it is essential 
we study, term by term, the analytic and asymptotic properties of the action as a function 
of these variables. It is easy to see that the action is an analytic function for all values of 
(pi. However, J as a function of ^ has several square root branch points, and logarithmic 
branching points and infinities. 

Note that a term arccos-u(2;) has branch points at u{z) = +1, —1, and at u{z) = oo. The 
associated branch cuts are usually taken to be (— oo, —1] U [1, +oo) . Since 



i\ogiu{z) + y-u(x 



arccos -uf^;) = — iloe;! u(z) + x/uiz)"^ — 1 



then we see that there are logarithmic singularities when u{z) = oo. The table below 
shows the logarithmic branching points and infinities associated with the dihedral angles. 
Furthermore, associated with the square root terms we have branch points a,t ^ = ^, C, = ^, 
and ^ = ^4- 
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Implementing the necessary branch cuts it is easy to conclude that the Riemann surface 
of the action is composed of an infinite number of sheets, each with a branch cut (— cx3,,^4]. 
However, as in [0], only three sheets will be relevant for the computation of the path integral. 
We define the first sheet Ci of I[C,] as the sheet where the terms in arccos(;z) assume their 
principal values. So the action in the first sheet will be formally equal to the original 
expression (|3.26|) . Note that with the first sheet defined in this way, for real ^ > ^4 the 
volumes and deficit angles are all real, leading to a real Euclidean action for ^ G [^4, +00) 
on the first sheet. On the other hand, when ^ is real and less than min{l/A, (3/4) in the 
first sheet , the volumes become pure imaginary and the Euclidean action becomes pure 
imaginary. For all other points of this first sheet the action is fully complex. 

Note that since we are only interested in geometries in which the boundary three-metric is 
positive definite, we must require that the volume of the boundary three-simplices be positive, 
which is equivalent to requiring ^~2 < P < ^^^ • Furthermore, since the simplicial metric 
in each four-simplex is real if and only if S, is real, then the simplicial geometries built out 
of these 4— simplices will be real only when ^ is real. Finally, computing the eigenvalues 
(Ai, A2, A3, A4) of the simplicial metric, gij^ we obtain 



Ai 
A2 



Se - (/3 + 1) - ^n-2{(3+l)f + {(3-lf 



Se - (/3 + 1) + ^n - 2(/9 + 1)]' + (/? - 1)^ 



A4 



^[l + /5-V5(/3-l)] 



(3.28) 
(3.29) 
(3.30) 
(3.31) 



It is easy to see that since ^ ^^ < (3 < "^^^ 1 then A2, A3, A4, are all positive and so the 
signature is Euclidean (+ + ++), when Ai > 0, i.e., when .^ > ^4 and Lorentzian ( — h ++), 



when Ai < 0, i.e., when ^ < ^4. So we see that for real ^ > .^4 we have real Euclidean 
signature geometries, with real Euclidean action, and for real ^ < min{^, j), we have real 
Lorentzian signature geometries with pure imaginary Euclidean action. 

When the action is continued in C, once around all finite branch points C, = 1/4, [3/4, ^^a, Csb, 
and ^4, we reach what shall be called the second sheet . It is easy to conclude that the action 
in this second sheet is just the negative of the action in the first sheet. 



I'[^,P,SA^ 



/^^[e,/3,^,0„0,] 



(3.32) 



Once in the second sheet, if we encircle the branch points in such a way that we cross the 
branch cut, (—00,^4], between 1/4 and [3/4, we arrive at what we shall call the third sheet. 



By doing this the terms arccos/ii, and arccos/i2, change signs, but J^ — 1/4 does not and 
so the action in this third sheet is different from in the first sheet 



Iiii[^,S,p,(l)i,(l)b] 



+ 



S \ ( 5 I 5 I 

— j|-^4/3 - /52(vr - 2arccos/i6) + -\J4^ - 1{-k - 2arccos/i5 



5 V 'C ~ l/4(27r + arccos hi + 2 arccos /12) 



5J pJ ^ ^"^^ ~ s-^ccos /i3 — 2 arccos h/^) 



4g >fi + 1^-/3^ + 3/3-1^^^515(0. - 0.)^/^ 



S 



h) (-)V/5 



iJ-p^ + W-i^i-U 



^^^,{^)m^^m 



(3.33) 



3.1 Asymptotic Behaviour of the Action 



In any discussion of the convergence of an integral over an infinite contour the asymptotic 
behaviour of the integrand, when | ,^ |^ 00, is essential. 
In the first sheet when ,^ ^ 00 the action behaves like 



l'[i,[3,S,<t>,,M ~ -v'(/5 + l)(-/92 + 3/5-1) 

1 m2/2 
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(3.34) 



where 






+ 



'(/5 + l)(-/52 + 3/5-ll 



(277 — arccos /i^ — 2 arccos /i^^ 



and 



arccos h'^ — 2 arccos h'^) 



i + ^^(W^ + # 



,2/2^ 



1-1 



(3.35) 
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ht 



/3(4 - /3) 

3/3-2 
^/?(4/3-l)(4-/3) 

2/3^ -2/5 + 1 



(3.36) 

(3.37) 
(3.38) 

v(4-/3)(4^-l)- -- (^-^^^ 

For the second sheet the asymptotic behaviour of the action is just the negative of that 
in the first sheet. For the third sheet we have 



4/3-1 



hf 



P 



■(3-2/?) 



rlll 



[tf3,SA^Ab] ~ y{f3 + l){-f3^ + 3(3-1) 



where 



1 m'^P 



—^(6 + b^,^.,n)\J^ 



(3.40) 
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((3 + l)(-(3^ + 3(3-l) 
J (3 {271 — arccos K^ — 2 arccos h' 



(27r + arccos h'^ 



2 arccos /i^) 



2/2 



1 + 
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2,2 I i2r2> 



■(40^/^ + 0,^/ 



-1 



(3.41) 



4 Classical Solutions 

The classical siinplicial geometries are the extrema of the Regge action we obtained above. 
In our minisuperspace model there are two degrees of freedom ^, 0j. So the Regge equations 
of motion will be: 

dl 

- = (4.42) 



and 

dl 



(4.43) 



They are to be solved for the values of C,, (pi, subject to the fixed boundary data S, (3, (p^. 
The classical solutions will thus be of the form ^{S,P,(f)h), and (p^i^S, P,(j)b)- The solution 
^{S, P, (ph) completely determines the simplicial geometry. 

Note that we shall be working on the first sheet. Of course, since on the second sheet 
the action is just the negative of this, the equations of motion are the same. And obviously 
every classical solution ^j{S,l3,4>b) located on the first sheet will have a counterpart C,jj of 
the same numerical value, but located on the second sheet, and so with an action of opposite 
sign, I[C,j{S, P,(pb)] = —I[C,jj{S, P,(pb)]- So the classical solutions occur in pairs. 

Inserting the second equation, ( |4.43|) , 



in the first equation ([4.42|) , we obtain a very long cubic equation on S for each value of ^, 
given fixed /3 and </>{,. This equation can then be solved numerically for ^, and by inverting 
the resulting solutions we obtain three branches of solutions ^ = C,d{S, l3,(pb)- For obvious 
physical reasons we shall accept only solutions with real positive 5*. In figure 2 we show one 
such solution for a fixed value of (3 and (pb- 

The general form of this solution is similar to the ones associated with other values of f3 
and (pb- We see that in general as ^ ^ ,^o = wim(l/4, /9/4), the value of 5* diverges to +oo. It 
is this branch that will represent the classical solutions for the late Universe. The fact that 
these solutions are Lorentzian means that semiclassically the wave function of the Universe 
should be dominated by the contribution coming from classical Lorentzian spacetimes like 
our own, as desired. 



As in previous cases the domain of solutions is divided by the hne S = Slj.^^{cl)^ where 



slM) = slM = <Pf) 



V'(/5 + l)(-/52 + 3/5-1) 
2 arccos K^) + J P{2ir — arccos h'^ — 2 arccos h 



{2tt — arccos h'^ 



1 + ^^K^l 



-1 



rriH^ 



where K — -^^p-. 

We see that for < S* < S^j.^^{cl), we will have: 

• Two pairs of real Lorentzian signature solutions C,j {S, P, (ph) = 0/ i^j Py 'Ph) ^ (— oo, .^o]; 
and S,i {S,(3,(f)b) = ^jj{S, f3,(f)b) E (— oo,.^o] with pure imaginary Euclidean actions. 

rp rp 

• One pair of real Euchdean signature solutions C,j {S,j3,(f)b) = C,ji{S, f3,(f)b) G [.^4, +oo), 
with real Euclidean action. 

For S > Sl^^^{cl) we have: 

• Only one pair of real solutions C,j{S, P,(f)b) = ^ii{S, ^.(pb) G (— oo,.^o] that correspond 
to Lorentzian signature simplicial metrics, and whose Euclidean actions, though sym- 
metric, are both pure imaginary. 

I[US,(3Ab)] = -mn{S,[3,<|>b)]=^I[US,C^<l>l>)] 

If we increase the value of m ot (pb, the value of Scrit decreases to zero and eventually the 
branch associated with the Euclidean regime vanishes. Furthermore we can see that S^j.^^{cl) 
as a function of the anisotropy j3, becomes infinite at the points of maximum anisotropy. 
See figure 3. 

5 Steepest Descent Contour 

After studying the analytical and asymptotic properties of the action we can now focus on 
the Euclidean path integral that yields the wave function of the Universe. 

^[S,f3,(f)b] = f De^0,e-^[^'^'^''^-'^''] (5.45) 



In our simplified models the result obtained from a contour C is not very sensitive to the 
choice of measure if we stick to the usual measures, i.e., polynomials of the squared edge 
lengths. So we shall take 

r/s- 9 

DiDct>, = — ^rf0, = ---didcP, (5.46) 

2'KIL'^ 2'KlH'^ 

As we have mentioned above there is as yet no universally accepted prescription for the 
integration contour to use in quantum cosmology. Following Hartle |0|, we shall accept 
that the main criteria any contour should satisfy are that it should lead to a convergent path 
integral and to a wave function predicting classical Lorentzian spacetime in the late Universe. 
The steepest descents contour over complex metrics seems to be the leading candidate. In 
the simplicial framework, complex metrics arise from complex-valued squared edge lengths, 
( ^J1) . The boundary squared edge lengths, S, and (3 have to be real and positive for obvious 
physical reasons. But the interior squared edge length, ^, can be allowed to take complex 
values. 

In general, a SD contour associated with an extremum ends up either at oo, at a singular 
point of the integrand, or at another extremum with the same value of Im{I). We have seen 
that when S is big enough the only classical solutions are a pair of real Lorentzian solutions 
iUS,PAb)Ai{S,P,4>b)) and (^jj{S,P,^b)MS,PAb)), where ^z = 0/ < mm(l/4,/5/4). 
They are located on the first and second sheets respectively, and so have pure imaginary 
actions of opposite sign. Given that their actions are different valued no single SD path can 
go directly from one to the other extremum. On the other hand given that 

m = [nrw 

and 

mi] = -mil] 

where * denotes complex conjugation, we see that the SD path that passes through ^jj 
will be the complex conjugate of the SD path that passes through ^j-. So the total SD 
contour will always be composed of two complex conjugate sections, each passing through 
one extremum, and this together with the real analyticity of the action guarantees that the 
resulting wavefunction is real. 

The SD contour passing through the classical solution {^c/(5', /3, 0;,), 0f (S', /5, 0;,)} is 



CsD{S,P,h) = |(e e R,<l>^) ■■ /m[/(5,e, A0„0,)] = /[ed(5,/3,0fe),0f(>5,/3,0,)]j (5.47) 

where R is the Riemann sheet of the action, and I{C,) = iliC,)- 

In figure 4 we show the result of a numerical computation of this contour for m = 1, 
0ft = 1, /3 = 1.5 and S = 50. 

The behaviour is similar to that for other values of the above variables. Going upward 
from the extremum, the SD contour proceeds to infinity in the first quadrant along the curve 
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The convergence of the integral along this part of the contour for any polynomial measure 
is guaranteed by the asymptotic behaviour of the action on the first sheet 
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As we move downward from the classical solution we immediately cross the branch cut 
and the contour enters the second sheet. The contour then proceeds to cross the branch cut 
once more, this time between 1/4 and /3/4, emerging onto the third sheet where it finally 
proceeds to infinity inside the first quadrant along the curve 
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As on the first sheet the convergence of the path integral along this section of the contour 
is guaranteed by the asymptotic behaviour of the real part of the Euchdean action along it 
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5.1 Semiclassical Approximation 

From our study of the SD contour for 5* > Scrit-, we can conclude that the range of integration 
in the neighbourhoods of the classical Lorentzian solutions give the dominant contribution to 
the integral, since the contribution of the other critical points, i.e., the infinities, is negligible 
given the asymptotic behaviour of the action. So in order to obtain the relevant information 
about the wave function of the Universe it is not really necessary to do the full computation 
of the integral along the SD contour. A semiclassical approximation based on the classical 
solutions found above will suffice. 

Since for S > Scrit these are real Lorentzian solutions with purely imaginary actions 
Ik = il[^ki^y Py4>b),(pi\S, (3,(f)b)], then for S > Scru, the semiclassical approximation is 
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This semiclassical approximation is specially good when the integrand is sharply peaked 
about the classical solutions (extrema), which is particularly true when the argument of the 
exponential is large. This will be the case for the large S of the late Universe and for the 
whole range of S when H^ = A/^/3 is sufficiently small as it is the case of our late Universe. 

In figure 5, we show the result of the numerical computation of one such semiclassical 
approximation. It is clear that the wave function exhibits the oscillatory behaviour that 
characterises the prediction of classical Lorentzian spacetime in the late Universe as desired. 

For 5* < Scrit the situation is not so simple as there are two pairs of classical Lorentzian 
solutions as well as the usual pair of Euclidean solutions. A semiclassical approximation 
could conceivably be based on any of these pairs. Furthermore, as we increase the values of 
(pb and m, the value of Scrit sharply decreases and so we can envisage a situation where the 
range of existence of these solutions practically vanishes. 

Computing the wave funtion for S = Scrit, ( figure 6), we see that as in the pure gravity 
case , 0, the wave function peaks for universes with large anisotropy, and small values of 
the scalar field. If as in the continuum case we consider a scenario of quantum nucleation 
of the Lorentzian universe at S* = Scrit, this result seems to favour the universes with larger 
anisotropy. 



6 Conclusions 

We have found that the resuhs obtained by us in 0] can be extended to anisotropic models. In 
particular, following 0, we have considered the simplicial minisuperspace based on the cone 
over the simplest anisotropic triangulation of the three-sphere, coupled to a massive scalar 
field, 0fc- The anisotropy is reflected in the existence of two different kinds of boundary edge 
lengths. However, since we admit only one kind of internal edge length, the minisuperspace 
is still two-dimensional as in the isotropic models considered in ^. This means that the 
path integral will not involve any new variables. We have found that for the late Universe 
the only classical solutions are pairs of real Lorentzian spacetimes like our own. We showed 
not only that there is a steepest descents contour going through these classical solutions but 
also that it yields a convergent path integral. These SD contours are very similar to those 
obtained in [^ , but their behaviour is slightly more complex due to the existence of a larger 
number of branch points in the action function. As in the isotropic models, the semiclassical 
approximation is quite good, specially when S is large or when H is sufficiently small, as it is 
in our late Universe. The computation of the semiclassical wavefunction shows its oscillatory 
behaviour, characteristic of the prediction of classical Lorentzian spacetimes, as desired. 

Combined with the results in and , the results of this paper show not only the ver- 
satility of the simplicial approach to minisuperspace models, but also that their predictions 
are generically in agreement with the results from similar continuum models. In all of them 
we find that by making the natural choice of the integration contour as the SD contour, the 
wave functions all predict classical Lorentzian spacetimes for the late Universe. 
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Figure 1: Anisotropic triangulation of the three-sphere 
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Figure 2: Classical solutions for /3 = 2.5 with a scalar field of mass m = 1 
and d)h = 1 




Figure 3: Variation of Scrit{d) with /3 e [^^, ^^ 




Figure 4: Steepest descents contour for (3 = l.b, H = 7 and scalar field of 
mass m = 1 when (pi, = 1 and S = 50. Coming from infinity in the first 
sheet it goes through the real extremum in the first sheet at .^ = 0.1337. 
Then proceeds downward onto the second sheet where it encircles the branch 
point 1/4, crossing the branch cut again, now between 1/4 and | = 0.375, 
at ^ = 0.3494. Thus entering the third sheet where it proceeds to infinity 
within the first quadrant. 
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Figure 5: Semiclassical wavefunction for S > Sent = 53.51 with scalar field 
of mass m = 1 and ^^ = 1, for (3 = 2.5, and H = 7. 




Figure 6: Wave function for S = ScritiP, 4>) with scalar field of mass m 
and H = 7. Note that the wave funtion peaks for large anisotropy 
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